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Abstract 

Starting from the divergence pattern of perturbation expansions in Quantum Field Theory and 
the (assumed) asymptotic character of the series, we address the problem of ambiguity of a function 
determined by the perturbation expansion. We consider functions represented by an integral of the 
Laplace-Borel type along a general contour in the Borel complex plane. Proving a modified form of 
the Watson lemma, we obtain a large class of functions having the same asymptotic perturbation 
expansion. Some remarks on perturbative QCD are made, using the particular case of the Adler 
function. 
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I. INTRODUCTION 



It has been known for a long time that perturbation expansions in QED and QCD are, 
under plausible assumptions, divergent series. This resuh obtained by Freeman Dyson for 
QED [1] was a surprise in 1952 and set a challenge for a radical reformulation of perturbation 
;heory. Dyson's argument has been repeatedly critically revised and reformulated since 
6] (for a review see also 3]), with the conclusion that perturbation series appear to be 
divergent in many physically interesting situations. To give the divergent series precise 
meaning, Dyson proposed to interpret it as asymptotic to F{z), the function searched for: 

oo 

F{z) ~ zeS, z^O, (1) 

n=0 

where iS is a region having the origin as an accumulation point, z being the perturbation 
parameter. Dyson's assumption of asymptoticity has been widely adopted. 

By this, the philosophy of perturbation theory changed radically. Perturbation theory 
yields, at least in principle, the values of all the F„ coefficients. This can tell us whether the 
series is convergent or not, but what we want to know is under what conditions F{z) can be 
determined from ([1]). If the series in ([1]) were convergent and the sign ~ were replaced by 
equality, the knowledge of all the F„ would uniquely determine F{z). On the other hand, 
there are infinitely many functions having the same asymptotic expansion (flj). 

This situation raises the problem of finding the "correct" or "physical" function F{z), 
using the knowledge of all the F^ coefficients (or, in a more realistic situation, the knowledge 
of several first terms only) of the series. The infinite ambiguity of the solution of this problem 
may be reduced if a specific field theory or model is considered allowing one to exploit some 
additional inputs of the specific theory. Useful information can be found in the papers Q], 
qI and the references therein. 

The objective of the present paper is to discuss the ambiguities of perturbation theory 
stemming from the assumed asymptotic character of the series. A class C of functions 
admitting a given asymptotic expansion is specified by the lemma of Watson, which we 
recall in section IIIII Watson's lemma, on the other hand, does not imply that C is the 
maximal class of that kind. In section HVl we present, and in section |V] we prove, a modified 
form of Watson's lemma. The modified lemma, which we refer to as Lemma 2 in this paper, 
allows us to show that the class C of functions possesing one given asymptotic expansion can 
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be, under plausible conditions, much larger than C. A discussion of Lemma 2 and its proof 
is placed in section jyll We discuss some applications in section IVIIt using as an example 
the Adler function 10| in QCD. 



II. PERTURBATION THEORY AND ASYMPTOTIC SERIES 

A. Pert ur bat ive approach 

A typical difficulty in physics is lack of exact solutions. To find an approximation, one can 
neglect some effects, which can then be reintroduced as series in powers of some correction 
parameter, z, written generically as ([1]), where F{z) is the function searched for. It is 
assumed that the expansion coefficients F„ are calculable from the theory. In most cases, 
however, only a few terms have been calculated and, in QCD, we seem to be near the limit 
of what can be calculated within the available analytical and numerical tools. 

A typical question in the 1950 's was whether a perturbation series of the type ([T]) was 
convergent or not. In many field theories and models, the large-order behavior of some 
subclasses of Feynman diagrams shows that the series is divergent, the coefficients Fn growing 
as n\ . But a sum can, under certain conditions, be assigned even to a divergent series. 
So, the crucial problem is: does ([T]) determine F{z) uniquely, or not? The answer depends 
on additional inputs and, also, on how the symbol ~ in ([T]) is interpreted. 

B. Basic properties of asymptotic series 

Definition: Let 5 be a region or point set containing the origin or at least having it as 
an accumulation point. The power series Yl'^=o Fnz"' is said to be asymptotic to the function 
F{z) as z ^ on 5, and we write Eq. ([1]), if the set of functions Rn{z), 



N 




(2) 



n=0 



satisfies the condition 




(3) 



for all N = 0,1,2,..., z 



and z G S. 



3 



We stress that an asymptotic series is defined by a different limiting procedure than 
the Taylor one: taking N fixed, one observes how Rn{z) behaves for z 0, 2; G 5, the 
procedure being repeated for all > integers. In a Taylor series, however, z is fixed and 
one observes how the sums X]^=o-^"'^" behave for 00. Convergence, a property of 

the expansion coefficients F„, may be provable without knowing F{z), to which the series 
converges. However, asymptoticity can be tested only if one knows both the Fn and F{z). 

The function F{z) may be singular at 2; = 0. The coefficients F„ in ([1]) can be defined by 

n-l 



^ri = lim — 



F{z)-Y,Fkz' 



(4) 



fc=0 

This definition makes sense whenever the asymptotic expansion ([T]) exists. To define F^, we 
do without the n-th derivative of F{z), z E S, which may not exist. 

Relation ([1]) does not determine F{z) uniquely; there may be many different functions 
with the same F„ coefficients. Note that the series ([1]) with all the Fn vanishing, Fn = 0, 
is asymptotic to many functions that are different from the identical zero. Let us denote 
the generic function of this type by H{z)] one example is H{z) = he~'^^^ with h ^ and 
c > 0. The expansion with all coefficients vanishing is asymptotic to he~'^^^ in the angle 
I argz| < 7r/2 — e, where e > 0. Then, F{z) and F{z) + H{z) have the same asymptotic 
expansion in the intersection of the two angles, in which the expansions of F[z) and H{z) 
hold. 

The ambiguity of a function given by its asymptotic series is illustrated in a more general 
formulation by the Watson lemma. 

III. WATSON LEMMA 

Consider the following integral 

<fo,c(A)= r e-^^"x^''^f{x)dx, (5) 
^0 

where < c < 00 and a > 0, /3 > 0. Let f{x) e C°°[0,c] and /*^^^(0) defined as 
limj.^o+ f^'^K^)- Let e be any number from the interval < £ < 7r/2. 

Lemma 1 (G.N. Watson): If the above conditions are fulfilled, the asymptotic expansion 

^_^^fk + (3\f('\0) 

k=0 



$o.(A)~if:A-^rf^ 

I— n \ 



a ) k\ 



(6) 



holds for A — s> oo, A G S's, where Se is the angle 

|argA|<|-£. (7) 
The expansion ^ can be differentiated with respect to A any number of times. 



For the proof see for instance 



ll|-|l3|. 



Remark 1: The perturbation expansion in powers of z discussed in the previous section 
is obtained by setting F{z) and I/2; in the place of $o,c('^) and A respectively. The formulae 
for F{z) corresponding to (JSj), ([6]) and (JTj) can be easily found. 

Remark 2: The angle ([7]) does not depend on a, (3 or c. 

Remark 3: The factor T^^^^ makes the expansion coefficients in ^ grow faster with 
k than those of the power series of f{x). 

Remark 4: The expansion coefficients in ([6]) are independent of c. This illustrates the 
impossibility of determining a function from its asymptotic expansion, as discussed in the 
previous section: the same series is obtained for all the integrals along the real axis, having 
any positive number c as the upper limit of integration. 

Below we shall display yet another facet of the above ambiguity, showing that under 
plausible assumptions the integration contour in the Laplace-Borel transform can be taken 
arbitrary in the complex plane. 

IV. A MODIFIED WATSON LEMMA 

Let G{r) be a continuous complex function of the form G{r) = r exp{ig{r)), where g{r) 
is a real-valued function given on < r < c, with < c < 00. Assume that the derivative 
G'{r) is continuous on the interval < r < c and a constant tq > exists such that 

\G'{r)\ < Kir^\ Tq < r < c, (8) 

for a nonnegative Ki and a real 71. 

Let the constants a > and /3 > be given and assume that the quantities 

A = inf ag{r), B = sup ag{r) (9) 

ro<r<c ro<r<c 

satisfy the inequality 

B-A<n-2€, (10) 
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where e > 0. 

Let the function f{u) be defined along the curve u = G{r) and on the disc |m| < p, where 
p > vq. Assume f{u) to be holomorphic on the disc and measurable on the curve. Assume 
that 

|/(G(r))| <f^2r^^ ro<r<c, (11) 

hold for a nonnegative K2 and a real 72. 

Define the function $[^'*(A) for < 6 < c by^ 

= r e-'^^^^-^r^G{r)f-'f{G{r))dG{r). (12) 



J r=b 

Lemma 2: // the above assumptions are fulfilled, then the asymptotic expansion 

k=0 

holds for A — > 00, A G 7^, where 



a ^ — ^ \ 



a J k\ 



V. PROOF OF LEMMA 2 



A. Proof 



(13) 



r, = {A : A = |A| exp(i(^), -'^-A + e<(p<'^-B-e}. (14) 



The conditions stated in section HVl assume implicitly that c > tq. We write: 

C(A) = Ci(A) + $S(A), (15) 
and define the new function G{r) by 

G{r) = — G(ro), for < r < tq; G{r) = G{r), for r > tq. (16) 
Since f{u) is holomorphic on the disc \u\ < p, p > tq, Cauchy theorem allows us to write 

<i(A) = <i(A), (17) 



^ This integral exists since we assume that f{u) is measurable along the curve u — G{r) and bounded by 
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i. e. the integral along the curved path can be replaced by an integral along the straight line 
u = r exp{ig{rQ)). Futhermore, on the disc |n| < tq, the function f{u) can be expressed in 
the form 

= E ^-ip-""' + ^ CM''^'- (18) 

k=0 

Then $o^q(A) can be written as 

^tlW = E 4o(^)^^Tr^ + r (-Ar<^e'"^(^°)) e^^(^°)(re*^('-«))^"V^(G(r))dr, (19) 
fc=i ■ 

where we defined 

I^^{X)= [ (re*^'("o))^-i+^exp(-Ar"e^"3(''°))e^3("«)rfr (20) 

for < 6 < c. 

It is useful to write 

4oW = CW-4oo(A), (21) 
since the first term, Iq^{X), can be trivially computed. We have 

Io,ooW = r^+'=-^exp(-ArV"f(""V^- (22) 



From the condition (ITOl) and the definition iHM it follows that, for A G 7^, one has 
Re [Ae^°5(^o)] > g. Therefore we can use the well known result 



x'-' exp(-/.x°) dx = — ^ r - , (23) 



which holds for Re/i > 0. Setting S = P + k and yU = A 6*"^'-'^°^ we obtain from ([22 



rfc 1 x-iM ^ (k + (3 



C(A) = -A-^r^ . (24) 



By inserting this expression in f[T^ and using f[T7|) . we write f[T^ in the form 



*g' (A) - i 5: A-^ r (^) ^ 4S(A) - E 4,.(A) 

A:=0 ^ ^ ■ fe=l 



/W(o) 



A;! 

+ / (re*^'("oY"^^^(G'(r))exp(-Ar"e*"^("o))e^^'("o)rfr. (25) 



We now proceed to the estimation of the terms in the right hand side of this relation. For 
^r^}c{X) we use the definition (fT2!) and notice that for A G 7^ the inequality 

Re(ArV"^('')) > |A|r"sin£ (26) 
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holds. Consequently, we have 
Using also ([8]) and (fTTI) . we obtain 

/•oo 

|$S(A)| < Kii^a / x^-^+^^+^^e-l^'""^'""da;, (28) 

■/ rn 



'ro 

which, with the transformation = t, becomes 

oo 







There exists such that 

|3+71+72~" 



< if^e*"'"" (30) 



on the interval [ro, oo). So, Eq. fl2^ leads to 

' '•"''^^ ^' - a(|A|-l)sin£ ■ ^ ^ 

The integrals /ro,oo(A) appearing on the right hand side of (|25l) can be estimated in the 
same way. By comparing ( l20l) with ( |T2i) . it follows that /^^ oo(A) are obtained from $r^L(A) 
by the replacements G'(r) re^^^'^°\ (3 ^ (3 + k and f{G{r)) 1. Setting (3 ^ (3 + k, 
Ki = K2 = 1 and 71 = 72 = 0, we obtain the bound 

poo 

l4loo(A)|</ x^-^+VI^I^^^'^^'fix, (32) 
which can be written, using (!29l) . as 

' ^' - a(|A|-l)sin£ ' ^ 

for a nonnegative, fc-dependent, constant 

Finally, for the last term in the r.h.s. of fl25|) we use the bound on rjv given in f|T8l) and 
obtain, by same procedure, the upper bound 

pro 

Cn / x^-^x^+^e-l^l^'^^^^^tix. (34) 
Jo 

This integral can be bounded using (123|1 with 5 = /5 + + 1 and /i = A sine, which leads to 

x^-^x^'+^e-^^^-^^^'dx = O ((|A| sine)-"^) . (35) 
The estimates ( ISTj) . (!33l) and (!35!l . inserted in the r.h.s. of (l25l) . show that 



/ 5 

N 



<.,.,-iE-^r(^)q^.o(,.r-^) (sa) 



fe=0 

for A ^ cxD, A G 7^. This completes the proof of ([T3|) . 
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B. Optimality 



There is a question whether the angle 7^ given in Lemma 2 can be enlarged. We show 
that the angle is maximal by proving that outside 7^ the relation fl5B]) does not hold. 

Let us take a special case g{r) = on [0,p], where < p < c < oo, g{c) = 7r/4. Let 
< g{r) < 7r/4 be fulfilled for r G [p, c], g{r) being a smooth function. We choose a special 
function / by taking f{u) = fi for \u\ < p and /(G(r)) = /2 for r G [p, c] where /i,/2 are 
two different, nonzero constants. We take a = (3 = 1. Certainly we have A = 0,5 = 7r/4. 
The assumptions of Lemma 2 are fulfilled and we obtain 

T, = {A : A = |A| exp^y;), -^ + e<^<^-£}, (37) 

where e is an arbitrary positive number. Define, as in the Lemma 2, 

<i(A) = r e-^^« /(G(r)) dG{r). (38) 

Now we choose a ray that lies outside 7^ , 

Li = {A : A = |A|exp(-i(7r/2 + (5))} (39) 

where < 5 < 7r/4. We shall show that the function (A) is unbounded along the ray 

This function can be written 

<i(A) = [/i + (/2-/i)exp(-Ap)]/A-/2exp[-Acexp(z(7(c))]/A. (40) 

Now for A G Li we have 

I exp(-Ap)/A| = exp(-|A|p| cos(-7r/2 - (5))/|A|, (41) 

which is divergent for |A| oo. The last term of fHUl) . 

|/2/A||exp[-Acexp(2(7(c))]| = I/2/AI exp[-|A|ccos(-7r/4 - 5)] (42) 

converges to zero as |A| ^00. It follows that Lemma 2 does not apply for A G Li, |A| ^ cxo 
(see (dl). 

Now we choose another ray, L2, which also lies outside %: 

L2 = {A : A = |A|exp(«(7r/4 + 5))}. (43) 
9 



Certainly 

I exp(-Ap)| = exp(-|A|pcos(7r/4 + 5)), (44) 
which converges to zero for A G L2, |A| cxo. For the last term we have 

|/2/A||exp[-Acexp(z(?(c))]| = I/2/AI exp[-|A| ccos(7r/2 + 5)], (45) 

which is divergent for A G L2, |A| ^ cxd. 

VI. REMARKS ON LEMMA 2 AND ITS PROOF 

Remark 5: Up to now we have assumed that c > p, where p is the radius of the 
convergence disc of f{u). If c < p, then the entire contour from to c can be deformed up 
to a straight line, and we have instead of fl25l) : 



c-.,-ii:-^r(^)/^^-i:.u. 

A;=0 ^ ^ k=l 



/W(0) 



+ I (re*^'("))'^-VAr(G(r))exp(-Ar°e*"f('=))e^^('=)dr. (46) 





The integrals I^ooW ^^"^ ^^^^ term in P6l) can be estimated as above, with the re- 
placement of To by c in (l33l) and These parameters do not appear in the asymptotic 
expansion in the l.h.s. of ( l36l) . which preserves its form. The difference is that the exponen- 
tial suppression of the remainder depends now only on c. The reason is that we can choose 
ro = c, and the conditions ([H]) - ( |TTi) are empty. 

Remark 6: Watson's lemma is obtained in the special case where the integration contour 
becomes a segment of the real positive semiaxis, i.e. g{r) = 0, and /(r) G C°°[0, c\. In fact it 
is enough to assume that /(r) G C°°[0, c'] with c' < c, where /(r) is measurable for r G (c', c) 
and bounded according to ( fTTl) . 

Remark 7: For applications in perturbation theory, we set A = in ( fT2l) . For sim- 
plicity, we take the particular case a = (3 = 1. Lemma 2 then implies that the function 

Ftc\^)= r e-^M/^/(G(r))rfG(r) (47) 

Jr=0 

has the asymptotic expansion 

00 

^o?(^)~E^'^'/^'nO) (48) 
10 
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FIG. 1: Three examples of integration contours in the u-plane. 



for z 



and z ^ Z^, where 




(49) 



Remark 8: The parameter e in the condition (ITO!) is hmited by < £ < 7r/2 — (5 — A)/2, 
but is otherwise rather arbitrary. Note however that the upper hmit of e may be considerably 
less than 7r/2, being dependent on the value oi B — A. This happens, in particular, if the 
integration contour is bent or meandering. We illustrate this in Fig. [Hby three contours C^, 
/c = 1,2,3 given by the parametrizations u = Gk{r) = rexp{igk{r)) with gi{r) = — 7r/30, 
g2{r) = 7r/2 — 7r/30, and gs^r) = 0.6 + 0.25r + 0.2sin97rr. We consider a = (3 = 1 and take 
the upper limit c = 1. Then we have B — A = {] for the constant functions gi and (72, where 
e can be any number in the range (0, 7r/2), while for g^ the difference B — A equals 0.59, 
and the condition ffTOj) gives < e < 1.27. As follows from f|T^ and can be seen in Fig. [2] 
left, if e is near zero, the domains 7^^^, k = 1,2,3, are large. However, Eqs. (ETi) and 
show that the bounds on the remainder of the truncated series are loose in this case. If on 
the other hand e is near its upper limit (see Fig. [2] right), the bounds are tight, but this has 
a cost in the fact that the 7^^^, k = 1,2, 3, the angles of validity of the expansion, are small. 

Remark 9: We note that the parametrization G{r) = rexp{ig{r)) does not include 
contours that cross a circle centred at r = either touching or doubly intersecting it, so 
that the derivative G'{r) does not exist or is not bounded. In particular, this parametrization 
does not include the contours 

(i) that, starting from the origin and reaching a value ri of r, return back to a certain 
value r2 < ri, closer to the origin, and 
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_] I I I I I } I I I I 1 u I I I I I I I I I I I I I I I I I I I 

-1 -0.5 0.5 1 -0.2 0.2 0.4 0.6 0.8 1 

Re X Re X 

FIG. 2: Three regions 7^^^ in the complex A plane given by Eq. (|14p . for the three contours 
shown in Fig. [TJ If e is small, 0.05 say (left panel), all the three angles are large, though 7^ 3 is 
smaller than the other two. For e = 1.2 (right panel), e is near the upper limit 1.27 imposed by 
condition (jlOp for A; = 3; as a consequence, the angle 7^,3 is very narrow. 



(ii) whose one or several parts coincide with a part of a circle centred at the origin. 

The contours under (i) can be included by representing the integral as a sum of several 
integrals along individual contours that separately fulfill the conditions of Lemma 2. In 
particular, if the integration contour returns back to the origin, the expansion coefficients all 
vanish. The cases under item (ii) can be treated by choosing an alternative parametrization 
of the curve. 

Remark 10: Extensions of the Watson lemma to the complex plane were considered 



briefly by H. Jeffreys in Ref. 



ll| . under somewhat different conditions. 



a) Identical is the condition requiring that the function f[u) is holomorphic on a disc, 

|m| < R. 

b) In Ref. HI], the integration is performed along a right angle in the complex w-plane, 
first along the real axis and then along a segment parallel to the imaginary axis, up to the 
upper point of the integration contour. So, it is assumed that the integral along this right 
angle exists, in particular, that the function f{u) is defined along it. In fact, / is assumed to 
be analytic in the region bounded by the initial contour and this right angle, except for some 
isolated singularities. The condition of Lemma 2 requires only that, outside the circle, f{u) 
is defined along the curve of the original contour of integration. These weaker assumptions 
imply that our approach is more general. 

c) Lemma 2 and its proof allow us to use the bounds (13T|) . (l33l) . and (135|1 to obtain an 
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estimate for the difference between (A) and its asymptotic expansion ([T3l) in terms of 
the constants introduced in the conditions (jH]), and (fTTj) . 

d) Finally, we remind the reader that the proof of Lemma 2 in subsection IV Al allowed 
us to obtain a remarkable correlation between the strength of the bounds on the remainder 
and the size of the angles where the asymptotic expansion is valid. Indeed, ( 13T|) . ( l33l) . and 
( 15^ depend on the parameter e, which determines the angles % and Z^, see ( IT^ and (H^ 
respectively. As was pointed out in Remark 8, the larger the angle of validity, the looser the 
bound, and vice versa. 



VII. REMARKS ON PERTURBATIVE QCD 

We take the Adler function [l^ 

V{s) = -s^ - 1 , (50) 
to discuss applications of Lemma 2, where n(s) is the polarization amplitude defined from 
z j d^xe*'^-^ < 0|T{V;(a;)K(0)^}|0 >= {q^q, - (?^.g')n(s). (51) 
Here s = q^ and is the vector current for light {u or d) quarks. 



In accordance with general principles 10|, ^^i^) is real analytic in the complex s-plane. 



except for a cut along the time-like axis produced by unitarity. In perturbative QCD, any 
finite-order aproximant has cuts along the time-like axis, while the renormalization-group 
improved expansion, 

V{s) = Di as{s)/TX + D2 {as{s)/Tif + D^, {as{s)/Txf + . . . , (52) 

has in addition an unphysical singularity due to the Landau pole in the running coupling 
g.J g). According to present knowledge, (|52|) is divergent, D„ growing as n\ at large n 



19|. 



A. Ambiguity of the perturbative QCD 

To discuss the implications of Lemma 2, we can define the Borel transform B{u) by 

B{u) = Y,b^u-^ bn = ^. (53) 

n>0 '^'^O ■ 
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It is usually assumed the series (153|) is convergent on a disc of nonvanishing radius (this 
result was rigorously proved by David et al. [20] for the scalar ip^ theory). This is exactly 
what is required in Lemma 2 for the Borel transform. 

If we adopt the assumption that the series fl5^ is asymptotic, Lemma 2 implies a large 
freedom in recovering the true function from its perturbative coefficients. Indeed, taking for 
simplicity a = f3 = 1 in (fT2|) . we infer that all the functions Vq^{s) of the form 



Po%(s) = ^ / B{G{r)) dG{r) , (54) 

PO Jr=0 

where a{s) = as{s)/TT, admit the asymptotic expansion 

oo 

Vl{s)r~.J2Dn{a{s)r, a,{s)^0, (55) 

n=l 

in a certain domain of the s-plane, which follows from (fT4|) and the expression of the running 
coupling a{s) given by the renormalization group. 

As mentioned above. Lemma 2 imposes weak conditions on B{u) and on the integration 
contour. Outside the convergence disc of ( l53i) . the form of B{u) (denoted f{u) in section 
IIVP is largely arbitrary, being restricted only by the rather weak conditions of Lemma 2. If 
the function B{u) defined by (!53|) admits an analytic continuation outside the disc (which is 
not necessary for Lemma 2 to apply), the analytic continuation can be used as input in the 
integral representation (!55l) . Then, if the contour passes through the analyticity domain, B 
say, more specific properties of Vq^{s) in the coupling plane can be derived, in analogy with 
the case of Borel summable functions (see 2ll|). 

The integral (15^ reveals the large ambiguity of the resummation procedures having the 
same asymptotic expansion in perturbative QCD: no particular function of the form J)q^[s) 
can be a priori preferred when looking for the true Adler function. 

The proof of Lemma 2 shows that the form and length of the contour, as well as the 
values of B{u) outside the convergence disc, do not affect the series (|55ll . contributing only 
to the exponentially suppressed remainder. As seen from the r.h.s. of ( l25l) or ( 146|) . the 
terms to be added to fl55|) are ^^^^(A) and Ir^^ooW^ where A = l/(/5oa(s)). The estimates in 
( |3T|) and (!33l) imply the remainder to (l55l) to have the form hexp{—d/ l3Qa{s)) ~ h (— A^/s)"^, 
where we used the running coupling to one loop. The quantities h and d > depend on the 
integration contour and on the values of B{u) outside the disc, which can be chosen rather 
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freely. So, contains arbitrary power terms, to be added to (l55l) .^ 



B. Optimal conformal mapping and analyticity 

In problems of divergence and ambiguity, the location of singularities of 'D[s) and B{u) 
in the a(s)-plane and, respectively, in the w-plane, is of importance. 

Some information about the singularities of B(u) is obtained from certain classes of 
Feynman diagrams, which can be explicitly summed and from general arguments 



based on renormalization theory 



15|, 



22l |. This analysis shows that B{u) has branch points 



along the rays u > 2 and u < —1 (IR and UV renormalons respectively). Other (though 
nonperturbative) singularities, for u > 4, are produced by instanton-antiinstanton pairs. 
Due to the singularities at u > 0, the series f l52|) is not Borel summable. Except the above 
mentioned singularities of B{u) on the real axis of the w-plane, however, no other singularities 
are known; it is usually assumed that, elsewhere, B{u) is holomorphic. 

To treat the analyticity properties of B{u), the method of optimal conformal mapping 



2J] is very useful. If the analyticity domain B is larger than the disc of convergence of (l53l) . 



one replaces (|53|) by the expansion 



B(u) 



E 

n>0 



Cr,/W 



where w = w{u) (with w{0) = 0) maps B (or a part of it) onto the disc 



f l56|) converges. fl56l) has better convergence properties than fl53|) : in 2J], Schwarz lemma 



(56) 



wl < 1, on which 



was used to prove that the larger the region mapped by w{u) onto \w\ < 1, the faster the 
convergence rate of ( l56l) . 



If w{u) maps t 
and, as shown in 



le whole B onto the disc \w\ < 1, fl56p converges on the whole region B 
24| . its convergence rate is the fastest.^ Then, the region of convergence 



of fl56|) coincides with the region B of analyticity of B{u). 



^ The connection between power corrections and Borel-Laplace integrals on a finite range was discussed 
also in 



23| 



^ This mapping is called optimal. In the particular case when B is the u-plane cut along the rays u < — 1 and 
u > 2, the optimal mapping reads w{u) = + u — — u/2)/{y/l + u+ -^/l — u/2) [24, 25|. Note that 
the expansion (|56p takes into account only the location of the singularities of B{u). Ways of accounting 



for their nature can be found in 



25 



26|. 



15 



Optimal conformal mapping allows one to express analyticity in terms of convergence. 
Inserting fl56|) into fl54|) we have 

T^Us) = ^ r e-W^ [w{G{r))rdG{r) . (57) 

This expression admits the same asymptotic expansion (155!) . However, containing powers 
of the variable w, it implements more information about the singularities of the true Borel 
transform B{u) than the series (!53l) in powers of u, even at finite orders. So, one expects 



that the finite-order approximants o 



physical function searched for |25l. l27|. 



( 1571) will provide a more accurate description of the 



C. Piecewise analytic summation 



The choice of the integration contour may have a fateful impact on analyticity. In 28|, 129] , 



two different contours in the w-plane are chosen for the summation of some class of diagrams 



(the so-called renormalon chains [16|, 113] ): one contour, parallel and close to the positive 
semiaxis, is adopted for a{s) > 0, another one, parallel and close to the negative semiaxis, 
is taken when a{s) < 0. As was expected, and later proved in [sO^, analyticity is lost with 
this choice, the summation being only piecewise analytic in s. Although this summation 
represents only a part of the full correlator, it is preferable to approximate an analytic 
correlator by a function which is also analytic, since analyticity in the s-plane is needed for 
relating perturbative QCD in the Euclidean region to measurable quantities on the time-like 
axis. 

On the other hand, as shown in 



31 



32| . the Borel summation with the Principal Value 
(PV) prescription of the same class of diagrams admits an analytic continuation to the whole 
s-plane, being consistent with analyticity except for an unphysical cut along a segment of the 
space-like axis, related to the Landau pole. In this sense, PV is an appropriate prescription. 

VIII. CONCLUDING REMARKS 

The main result of our work is a modified form of Watson's lemma on asymptotic series. 
This result, referred to as Lemma 2, holds if the function f{u) (which corresponds to the 
Borel transform of a QCD correlator) is analytic in a disc and is defined along a contour. 
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on which it satisfies rather weak conditions, which are specified in section IIVI 

Our result emphasizes the great ambiguity of the summation prescriptions that are al- 
lowed if the perturbation expansion in QCD is regarded as asymptotic. The contour of the 
integral representing the QCD correlator and the function B{u) can be chosen very freely 
outside the convergence disc. 

We kept our discussion on a general level, bearing in mind that little is known, in a 
rigorous framework, about the analytic properties of the QCD correlators in the Borel plane. 
If some specific properties are known or assumed, the integral representations discussed in 
ths paper will have additional analytic properties, known to be shared also by the physical 
amplitudes. 

Lemma 2 proved in this paper may also be useful in other branches of physics where 
perturbation series are divergent. 
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